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A formalism to study light fermions at non-zero chemical potential in the instanton 
ON ! medium is presented. It is then applied to investigate chiral symmetry restoration and 
diquark condensation at finite baryon density. 
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o . 

OO ■ 1. INTRODUCTION 

ON ; 

^ ■ The fact that instanton-induced interactions produce attraction not only in the qq 

channel (leading to the spontaneous chiral symmetry breaking) but also in the qq channel 
(potentially leading to diquark condensation) has been first realized in ref . ]1| . It becomes 
especially clear when the number of colours is N c = 2, where instanton-induced interac- 
tions possess a global SU(4:) symmetry H (often referred to as Pauli-Giirsey symmetry). 
For this reason qq and qq condensates belong in fact to one phase: one condensate can be 
rotated to another along the Goldstone valley 0. 

The possibility of diquark condensation (in analogy to superconductivity) at any iV c , as 
induced by instantons, has been studied in ref. 0. Only a metastable diquark-condensed 
vacuum has been found at N c > 2 and chemical potential /i = 0. Furthermore, already at 
N c = 3 the scalar diquark appears to be unbound in the vacuum, indicating that our world 
is in a sense closer to the idealized N c — > oo than to the N c = 2 world. Parametrically, 
the diquark mass is ~ 1/p ~ 1 GeV where p is the average instanton size as explained 
below. It means that a scalar diquark correlation function should decay with the exponent 
corresponding to the 'constituent' quark threshold 2M(0) ~ 700 MeV, which seems to 
be supported by recent lattice measurements [Q. Nevertheless, it has been suggested in 
that qq condensates could be found as metastable states in heavy ion collisions and in 
astrophysics. 

Quite recently diquark condensation owing to instantons at fi ^ has been estimated 
by the Stony Brook and Princeton/MIT groups HH, see also these proceedings. Below 
we present a formalism to study it in a more systematic way. A theory of light quarks 
in the instanton vacuum, though not completely trivial, is rather well developed, and 
it is straightforward to generalize it to /i ^ 0. We will avoid some of the unnecessary 
approximations made in ||§]. 
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2. QCD INSTANTON VACUUM 



The Euclidean QCD partition function for Nf fermions with a common chemical po- 
tential n is 



Z = j DA^exp {-J ^ J D^j Dtf exp £ J 4>}(if+ im f - ip^ipf 



It is understood that the theory is regularized in some way at large momenta, for example, 
by the Pauli-Villars method. It is also convenient to normalize the partition function to 
the free one. The integral over fermion fields can be thus understood as 



Dip Dip^ exp 



Det 



n 



det (if + im f - ^74) det (2^ + iM) 
det(«y+ iM) det(i# + im f ) 



(2) 



where M is the Pauli-Villars regulator mass playing the role of the ultraviolet cutoff. 

The main hypothesis of the instanton vacuum is that the partition function (|IJ) is 
mainly saturated by the ensemble of instantons and antiinstantons (J's and J's for short), 
plus quantum fluctuations about them. Since the ensemble of J's and J's is not, strictly 
speaking, a solution of the Yang-Mills equation of motion (it is only in the infinitely dilute 
limit), one has to do something more clever than the saddle-point integration over A^. 
The most straightforward procedure is to use the Feynman variational principle in order 
to get the lower bound for the partition function (fj) 0. Writing the general Yang-Mills 
field as a superposition of N + J's and N_ J's plus a presumably weak quantum field B^, 



(3) 



where £ is a set of 4N C collective coordinates for each instanton (namely center, size and 
colour orientation), one can show that the partition function ([!]) can be approximated by 



^>E 



NJN. 



Jdi J(0 exp[-C/ int (e)].Det. 



(4) 



Here J is a Jacobian which is due to passing to the collective coordinates and Ui nt is 
the instanton interaction. Actually, the separation of the dependence on the collective 
coordinates £ into J and U int is ansatz-dependent since the precise meaning of the collec- 
tive coordinates depends on the functional form of the field. However, their combination 
should be ansatz-independent, at least if one calculates both quantities with infinite pre- 
cision. We shall assume the simplest sum ansatz 

It has been shown in ref. (see also a more recent paper ||) that the partition function 
(f|) (without the fermion determinant) leads to the stabilization of the grand canonical 
instanton ensemble. The basic characteristics of the medium is the average instanton size, 
p, and the average separation between near neighbours, R, or the instanton density at 
equilibrium, N/V = 1/(R) 4 . Since the transmutation of dimensions takes place in this 
approach, both quantities are expressed through the only dimensional parameter one has 
in the pure glue theory, namely Aqcd- The 2-loop calculations performed in J7||| give 
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Taking = 280 MeV one finds p « 0.35 fm, R « 0.95 fm, p/R « 1/3. This small 
ratio has been previously suggested on phenomenological grounds by Shuryak ||. The 
smallness of the pj R ratio implies that the packing fraction of instantons in the vacuum, 
i.e. the fraction of the 4d volume occupied by the balls of radius p, is quite small: 

7 2 P 1/ 2 i? 4 ~ 10' U 

2.1. Instanton Packing Fraction as Algebraic Parameter 

After the transmutation of dimensions takes place all physical quantities are propor- 
tional to Aqcd in appropriate powers with numerical coefficients of the order of unity, as 
it should be, generally speaking, in a strong interaction theory. Therefore, if one gets a 
numerically small (or large) dimensionless quantity it must be due to some exceptional 
reason. In our case the reason why one gets a numerically small packing fraction (|^) can 
be traced back to the 'accidentally' large value 11/3 of the coefficient in the asymptotic 
freedom law. The same analysis applied to the 2d CP N model also possessing instantons 
shows that there are no such numerical tricks there: instantons of that model seem to be 
strongly overlapping, and there appears no small packing fraction parameter. 

Though p/R ~ 1/3 is but numerically small, one can treat it as a formal algebraic 
parameter, and develop a perturbation theory in it. In principle, one can change the 
parameter p/R by taking more quark flavours, or taking nonzero temperature or including 
Higgs field which cuts off integrals over the instanton sizes. 

2.2. Separating High and Low Fermion Eigenmodes 

We now switch on Nf flavours of light quarks represented by the fermion determinant 
(|2T) in the backround field of the instanton ensemble. The crucial feature of fermions in 



the one-instanton background is that there is an exact zero mode ||10|1 . In the infinitely 
dilute limit of instantons that would lead to the nullification of the fermion determinant. 
Though the packing fraction of instantons is small it is, however, nonzero, and one can 
build a systematic theory of light quarks in the instanton vacuum by treating p/R as a 



formal small parameter |]TT]^T3 



An important step is to separate contributions of high- and low-frequency parts to the 
fermion determinant. The role of the high-frequency part is to renormalize somewhat the 
one-instanton weight or fugacity. The low-frequency part is physically much more inter- 
esting: it is responsible for chiral symmetry breaking in the vacuum and its restoration 
at high densities and/or temperatures. 

The separation is made in respect to an auxiliary mass Mi lying inside a parametrically 
wide interval, 1/R <C Mi <C 1/p. One rewrites the fermion determinant (@) as a product, 

detgy + im f - ipj 4 ) det (iff + iM) 
det(2y + iM) det{i(fi + im/j 

where 

_ det(?y+ iM x ) det(i$ + iM) _ det (if + imf- ipj±) det(i$ + iM x 

etHi9h ~ det(tf+tM) det^ + iM^ 6ti< ™ ~ det{i 1 f+ iM\) det(i$> + irrif) 

These ratios are arranged so that Detnigh gets a contribution from fermion modes with 
Dirac eigenvalues ranging from Mi to the Pauli-Villars mass M, while Det Low is saturated 
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by eigenvalues less than M\. The product of the two determinants is, by construction, 
independent of the separation scale M\. However, we are going to treat both of them 
approximately. The fact that the actual dependence of the product on M\ turns out to 



be extremely feeble []1T| in the wide range of M 1 serves as a check of the approximations. 

In considering Det High) one assumes that Mi is large enough so that Det^^ can be 
factorized into a product of determinants computed in the background of individual in- 
stantons. Corrections to the factorization are of the order of (MiR)~ 2 <C 1. 

In considering Deti ow , one assumes that Mi is small enough so that only the (diagonal- 
ized) would-be zero modes of individual instantons can be taken into account. Corrections 
from nonzero modes of individual instantons are of the order of (Mip) 2 <C 1. 

Taking nonzero chemical potential fi does not seriously change the above logic. At 
fip <C 1 (which will be our domain of interest since all important phenomena happen 
inside this range) the [i term should be included into the 'low' part as in eq. (0). 

2.3. Three Ways to See Chiral Symmetry Breaking 

In refs. [|i~T|-|13|j three seemingly different but actually equivalent ways of computing 
DetLo™ have been developed, thus leading to several languages in which chiral symmetry 
breaking by instantons can be understood and mathematically described. 

The first method is based upon diagonalizing the overlaps of the would-be zero fermion 
modes of individual instantons ITTIJT^l. One can write 




Det Lmi) = exp \ ^ I dX z/(A) In ^ 2 - ^ \ (8) 



where is the spectral density of the Dirac operator in the field of the instanton 
ensemble. It is obtained from diagonalizing a matrix made of the overlap integrals of the 
would-be zero modes, and then averaging over the ensemble. The nonzero z/(0) signals 
chiral symmetry breaking: it is actually due to the derealization of the zero modes owing 
to the "hopping" of quarks from one instanton to another. 

The second method is to find the quark propagator in the instanton ensemble. In 
the zero-mode approximation (justified when the packing fraction (|6]) is small) one can 
write a closed equation for the quark propagator averaged over the instanton ensemble 
0,14]. The propagator appears to be that of a massive fermion with a momentum- 



dependent mass M(p) related to the Fourier transform of the zero mode. Using the main 
characterictics of the instanton vacuum @ one gets M(0) ~ 350 MeV which is what one 



expects for the so-called constituent quark mass. Parametrically, M(0) ~ y N/VN c p, 
meaning that the constituent quark mass is small in the packing fraction of instantons. 

Integrating the quark propagator over the 4-momentum one gets the chiral condensate 
{ipip) ~ —(250 MeV) 3 being also close to its phenomeno logical value. 

The Detiow can be expressed through the propagator S(x,y) as 

f Ml 

Detz,o«j = exp / i dm' Tr [S(x, x, m!) — Sq(x, x, m')\ (9) 

J m 

where So is the free propagator. Comparing the integrand in eq. (D with that of eq. (|) 
one finds the Wigner semicircle spectral density of the Dirac operator averaged over the 
instanton ensemble ||12|| , which can also be obtained directly from the first method fl5 |. 
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Both methods can be directly generalized to nonzero chemical potential, since at /x 
there also exists an exact normalizable zero fermion mode M,[[7j • However, in this paper 



we use the third method suggested in |13| which is also easily translated to /x 7^ 0. 



3. EFFECTIVE FERMION ACTION AT /x ^ 

Rather than calculating fermionic quantities in the background of instantons and then 
averaging over these configurations, one can perform the averaging procedure first. This 
generates instanton-induced interactions in the form of a vertex involving 2Nf quarks. It 
is a direct consequence of the appearance of the quark zero modes, one for each flavour, 
and the resulting interactions are naturally of the 't Hooft form Jl0| . 

The guiding principle is to write a fermion partition function reproducing the prop- 
agator at low momenta as well as the fermion determinant. This derivation has been 
developed in detail in previous publications |L3| and its extension to finite chemical po- 
tential is straightforward. It should be stressed that correlations between instantons 
induced by fermions are inherent in this approach; as to correlations induced by gluons, 
they are effectively taken care of by the use of the variational principle [0] resulting in 
the effective size distribution. For simplicity we freeze all the sizes at the average value 
p, but average explicitly over randomly positioned and oriented instantons. 

The partition function which produces the necessary propagator is of the form 



Z = Det 



Low 



Dip Dtp^ exp 



f 



N f 



Y, 



N- 



N f 



V 



V 



(10) 



The pre-exponential factors contain the instanton-induced interactions between fermions, 
and for Nf flavours are the specific nonlocal 2A r ^-fermion vertices 



Here we use the notation x = x^a^, where the 2x2 matricies = (±z<7, 1) decompose 
the Dirac matricies into chiral components, and it is understood that /x written as a four- 
vector is ix a = (0, /x). Note that the zero mode in the field of an (anti) instanton couples 
to that of a (right-) left-handed quark. By we denote the conjugate zero mode. As 
noticed in [|16|], $/(x, /x) = $}(x, — /x), where the dagger means hermitean conjugate. 

Fermion operators in the pre-exponent are not convenient; these operators can be raised 
into the exponent with the help of a supplementary integration over a pair of Lagrange 
multipliers, denoted A±: 

Z = J d\ + d\_ J exp I J d 4 x - x/174)^ + X+Y+ f + \_Y Nf 



+N+ In 



N + 



1 + i\L In 



A3 7 



(12) 



Indeed, integrating over A± by the saddle-point method one recovers eq. fllO| ) (the saddle- 
point integration becomes exact in the thermodynamic limit N± —>■ 00). Through this 
procedure we obtain a purely exponential integrand which is the required effective action. 
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It is very important that the overall strength of the 2iV/-fermion interaction, whose role 
is played by \±, is not fixed once and forever: actually its value is found from minimizing 
the free energy after integration over fermions is performed. Therefore the strength of 
the interaction depends itself on the phase the fermion system assumes. As will be shown 
below, this is the mechanism by which chiral symmetry is restored at large chemical 
potential. 

Although one can investigate the phenomena of chiral symmetry breaking and its 
restoration for any number of flavours, we shall restrict ourselves to the case Nf = 2 
corresponding to a system of chiral up and down quarks. Furthermore, we will assume 
the CP invariant case of 9 = 0, which requires N + = i\L = N/2 and hence A + = A_ = A. 
For practical applications it is favourable to use the Fourier-transformed expressions of the 
quark zero modes, which are written explicitly in the Appendix. These complex functions 
of the four-momenta and chemical potential determine the (matrix) form factors attached 
to each fermion leg of the vertex: 

F(p,t*) = (p + w)~<p(p,v) + , F*(p, -fi) = <P*(p, -fi)~(p + in) + . (13) 
With such definitions, the interaction terms may be written in momentum space, 

KY = K J (2tt) 5 (pi +p 2 -h- k 2 ) 

■V'k« <a (ft)^(ft./')£^ fa ^^e na o 9 ^(*a. 2P2 (h) , (14) 

with a similar form for \ JY~ which carries right-handed quarks. The first indices on the 
fermion operators refer to flavour (1 or 2 explicitly), the Greek to colour (1 . . . N c ), and 
the last denote spin (1,2). This formulation of the effective interaction retains the full p 
and \i dependence of the zero modes, as opposed to the approximated treatments in other 
recent works [0,0. 

4. COMPETITION BETWEEN qq AND qq CHANNELS 

Since the instanton-induced interactions ( |14"D support both qq and qq condensation, it 
is necessary to consider the two competing channels simultaneously. This means that 
one must calculate both the normal (G) and anomalous (F) quark Green functions. A 
colour /flavour /spin ansatz compatible with the possibility of chiral and colour symmetry 
breaking is 

(¥ a \p)^UM) = Sf 9 ^ S i(P)) for «. P = 1. 2 > 
(V /ai (p)^(p)> = #5?&(P)$ for a, P > 2 , 
(i>l ai (p)i>f j (-p)) = (#(P)#'(-P)> = ef'fPW&Fip) , (15) 

where [7] refers to some generalized direction(s) in colour space, and it is this set of N c — 2 
indices which signals the breaking of colour symmetry. In the particular case of N c = 3, 
where the colour symmetry is broken as SU(3) — ► SU(2) x U(l), we will by convention 
take [7] = 3; for N c = 4 one can take [7] = 34 and so forth. In the event of colour 
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Figure 1. Gorkov equations to first order in A. 

symmetry breaking, the standard propagators (and ensuing condensates) will lose their 
colour degeneracy and the separation of S(p) into Si(p) and S 2 (p) becomes necessary; 
otherwise the Schwinger-Dyson equations do not close. 

Written in the chiral L, R basis, the 4x4 propagators S\^{p) are of the form: 



S(p) 



G(p)i z( P )s ( P y 

Z(p)S (p)- G(p)l 



So(p) ± = (p + ifi) 



(16) 



Here the off-diagonal, bare propagator is modified by the scalar functions Zi >2 (p), and is 
augmented on the diagonal by the scalar G\^{p) which if nonzero break chiral symmetry. 

Using the instanton-induced interaction ( |I~4"D one can build a systematic expansion 
for the F, G Green functions in the 1/N C and p/R parameters. In the leading order in 
both parameters we restrict ourselves to the one-loop approximation shown in Fig. I[ It 
corresponds to a set of self-consistent Schwinger-Dyson equations (called Gorkov equations 
in the case of superconductivity with F^O). An important /^-dependence enters through 
the form factors in eq. fll4|) . With F(p) = F*(p) these diagrams lead to the set of five 
algebraic equations for the scalar Green functions Z X;2 , G x>2 and F: 

Z x (p) = 1 - G x {p)A(p^)g x - F(jp)B(jp,rif 
Z 2 (p) = 1 - G 2 (p)A(p, [x)g 2 
Gl(p) = Zifytfafa fi)(p a (p, n)gi 
G 2 (p) = Z 2 (p)<p a (p, ti)ip a (p, fi)g2 

F (p) = Z(-p)ip a (p, n)ip a (-p, p)f . (17) 



The numerical factors g x , g 2 , and / will be defined below and the functions 



AfafJ.) = (p + i/x) a (P + ¥)a<Pf}(P, V)<Pp(p, y) , 

B{p,p) = (p 2 + /i 2 )^/3(p,/i)^/3(-p,/i) + (p + ifj,) a <p a (p,fj,)(p- ifj,)pipfi(-p,fj.) 

-(p + ifl)a<Pa(-p, y)(P ~ ifl)ptpp(p, H) , (18) 

are the form factors which arise from the zero modes (see the Appendix). At \x = we 
have A(p, 0) = B(p, 0), but for any finite [i the direction of the momentum flow through 
each vertex leg is critical. 

The condensates g±, g 2 , and / are the closed loops contributing to the quark self-energy. 
They are found by integrating the appropriate Green function, modified by the vertex form 
factors, over an independent momentum: 



8 



These appear in the particular combinations 



9i 
f 



5 - U + ( 2iV c - 5 



92 



, 92 = 2(2- j^jg 1 + 2(N c -2)g 2 



2 1 



1 



./ 



(20) 



Although the integrands above are complex, all imaginary parts are odd in p^ and thus 
vanish under integration. As the 91,2 are measures of chiral symmetry breaking, these 
act as an effective mass modifying the standard quark propagation. On the other hand 
the diquark loop / plays a different role, that of twice the single-quark energy gap (con- 
ventionally denoted A) formed around the Fermi surface. The Fermi momentum, in the 
absence of chiral symmetry breaking, will remain at pf — /i regardless of the magnitude 
of/. 

After determining the five scalar functions through solving eqs. QT7D and inserting these 



solutions into eqs. (|I~9"D we find the coupled equations for the condensates themselves: 



9i 
f 



A/ 



d 4 k a(l- (3f 2 + a*9l) 



d 4 k (3(1 - (3f 2 + agj) 



92 



A#2 



d 4 k 



a 



N*-lJ (2vr)4 (1 + agf)(l + a*gf) - (3 2 f 4 ' 
where yet another pair of functions has been introduced, 

a(p, n) = A(p, n)(p a (p, l*)<p a (p, n) , Pip, ft) = B(p, fi)(p a (p, fi)<p a (-p, fi) 



-2 ' 



(21) 



(22) 



Note that while j3 is real, the funtion a is complex. As usual for gap equations, there is 
the possibility of a solution where some (or even all) of the condensates vanish. 

The magnitudes of these condensates are not yet determined, as this requires fixing the 



coupling constant A. This is done through minimizing the partition function (12) and 
leads to 



N 
V 



X(Y + + Y-}. 



(23) 



On the left-hand side is the instanton density, which we will here take to be fixed at its 
vacuum value, although in principle will have some correction due to the finite quark 
density. Evaluating the right-hand side requires calculating the one-vertex contributions 
to the free energy, which in this case includes two-loop, figure-eight type diagrams formed 
by joining the four fermion legs into two pairs as shown in Fig. [| This can be written 
concisely in terms of the condensates: 



N 
V 



X(Y + + Y-) 



4(iV c 2 - 1) r 



A 



20i& + (N c - 2)g 2 g 2 + Iff 



(24) 



This equation for A and the definitions of the condensates ( pi|) comprise a system of 
equations which can be solved for all quantities. 

Once this has been achieved the chiral condensate may be calculated, which is a closed 
trace over the quark propagator. This is distinct from the condensates 9i 2, however it 




Figure 2. Contributing diagrams to Q at order A 1 . 

does contain the same properties under p 4 reflection which guarantee it is real. It is most 
simply written 

- (W)Mink = l^H)Eucl = *J-0yi^ S (p) = 4 / 7^4 l^l(p) + (N c - 2)G 2 ( P )\ .(25) 

As a consequence of having two different modes of quark condensation, one obtains mul- 
tiple solutions for the characteristics of the quark medium at any fixed chemical potential. 
Specifically, there will be a competition between the following phases: 

(0) No condensation of any type; g\ — g 2 — f — 0. 

(1) Chiral symmetry breaking without diquark condensation. This, the standard vac- 
uum scenario, has the signature g\ — g 2 ^ 0, / = 0. 

(2) Diquark condensation without chiral breaking, or / ^ 0, g\ = §2 = 0. 

(3) A mixed phase wherein both symmetries are broken, and g%, g 2 , and / are all distinct 
and finite. 

Phase (0) requires that A vanish, and such a solution was never found in our calcula- 
tions. To resolve between the remaining, symmetry-breaking phases, the free energy is 
minimized. Consistent with the evaluation of the Green functions, it is calculated to first 
order in A. Repeating the calculation of figure-eight diagrams and recalling the explicit 
dependence on A in eq. ( |T2"| ) we obtain the free energy 

1> [** + <* - + ^ .<*> 

Here V3 is the three-volume while V is the Euclidean four-volume, and Qq is the free 
energy for a gas of free quarks. The last two terms are precisely the quantity which must 
vanish under the saddle-point determination of A, and thus we have 

V 3 Vs V \N/V) k ' 

Thus the phase which features the lowest coupling A is the thermodynamically favoured. 

Numerically, no solutions were found which favour phase (3), the mixed case. Thus 
we concentrate on the distinction between phases (1) and (2), the former being marked 
by nonzero g = g\ = g 2 , the latter by nonzero /. Each case leads to its own value of A 
through solving eq. (^) for any given \i. At certain value of \i = fi c the solutions for A 
corresponding to the two phases cross; this is the point where the phase transition occurs. 
This point is defined by the condition 



?M = V 2 ' (28) 

When the ratio between the condensates is less than the constant on the right, phase (2) 
is favoured; otherwise it is chiral symmetry that is spontaneously broken in phase (1). 
Calculations were carried out for various iV c , taking the values for N/V and p specified 
in Section 2. 
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V (MeV) 

Figure 3. Condensates for N c — 3 as a function of fi. Shown are the effective quark mass 
M(0), the quark condensate (tptp) 1 ^ 3 , and the diquark energy gap per quark f/2. 

4.1. The Case of N c = 2 

In this case it is obvious that colour symmetry is not broken by diquark formation, 
which here correspond to colour-singlet 'baryons', and hence there is only one possible 
chiral condensate gi = g. Not so obvious is that at \i — the colour and flavour SU(2) 
groups are arranged into the higher SU(A) symmetry The instanton vacuum accounts 
for this symmetry 0, and in the context of the analysis here this corresponds to f 2 + g 2 
being the only discernable quantity in the gap equations (f2l"l). Numerically, we find 
Vf 2 + g 2 = 147 MeV. 

At finite /i, however, the SU(4) symmetry is explicity broken and one finds f < g. 
Since in this case the critical ratio of eq. ( |28| ) is unity, we conclude that for any finite 
density the N c = 2 world prefers diquark condensation to chiral symmetry breaking. This 
finding is in agreement with the reasoning of ref . || and the lattice results of ref . . 

4.2. The Case of N c = 3 

For three colours, the critical ratio of the condensates (|28|) is \f?>. At [i = we find f / g = 
165 MeV/65 MeV > v^- It means that at low values of /i the coupling constant A is smaller 
in the chiral broken phase, hence this phase is energetically preferred. Furthermore, the 
effective quark mass is the reasonable M(0) = g = 346 MeV. We thus recover the standard 
chiral symmetry breaking at low baryon density. 

With increasing chemical potential both / and g condensates in their respective phases 
are reduced, however as /i surpasses the effective quark mass M(0) around 300 MeV, the 
chiral-breaking g stabilizes while the colour-breaking / continues to decrease in phase (2). 
At /i ~ 340 MeV condition (|28|) is reached, and therefore for this and greater chemical 
potential chiral symmetry is restored and the 'colour superconducting' phase is realized. 
At the transition point, the quark gap is f/2 = 115 MeV and slowly decreasing with rising 
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density. Numerical results for the characteristic quantities of each phase are plotted in 
Fig. |3|. It is noteworthy that the chiral condensate, (ipip), is distinct from the effective 
quark mass, g, in that while the latter decreases with p the former remains practically at 
its vacuum value. For p> p c they both vanish. 

5. CONCLUSIONS 

We have formulated the effective low-energy fermion action induced by instantons at 
nonzero chemical potential. In the resulting fermion vertex we have retained the full 
dependence on both momentum and chemical potential, which arise from the would-be 
zero modes. In this respect we differ from other studies such as the random matrix model, 
NJL models, and alternative instanton approaches. 

The effective action leads to a competition between two phases, one of chiral symmetry 
breaking and another characterized by diquark condensation. It was studied by solving a 
coupled system of gap equations to first order in the instanton density. For two massless 
flavours and three colours, it was found that the usual broken chiral symmetry is restored 
through a first order phase transition, replaced by colour breaking due to the formation 
of a diquark condensate. This occurs at a critical chemical potential of ~ 340 MeV. 



APPENDIX. Fourier Transforms of Fermion Zero Modes 

The use of the exact fermion zero modes in the momentum space tremendously simplifies 
all calculations. The starting point is the exact fermion zero mode in the field of one 
(anti)instanton in x space [16,17|, which we cite for arbitrary instanton position z, size p 
and orientation given by rectangular N c x 2 matrix U: 



[$(x - z) 



IL(x) = 1 + 



„2 



X4-Z 4/ 



U(x-z) 



' e -A*(*4-**)A(a; - z,p)\ 1±75 



Tl(x - z) 



xj + r 2 + p 2 



A(x,/i) 



x\ + r 2 



X.\ 



cos(yur) H sm(pr) 



e> k U%.(29) 



(30) 



Its Fourier transform is defined as p) = i J d A x e~ ip x Q(x, p) and has the structure 
^(P)A 4 ) = la^Pa{p, p)- The Lorentz symmetry is broken at finite p, and the components 
become 



<^4(P4,p;/i) 



Tip 



y {(P - A* ~ m) [( 2 P4 + ifi)fl- +i(p~p- ZP4)/2-] 

+(p + v + m) [(2^4 + iy)h+ -i(p + p + ipa)M } 



T^p 2 Pi 



4 p 2 



(2p - p){p- p- ipA)h- + (2p + p){p + p + m)h- 



+ 
+ 



2(p - p)(p- p- ip A ) --(// + ip4)[pl + (p~ ^) 2 ] 

P 

2(p + A*)(p + P + ipi) + ~{p + iPi)[pl + (p + p) 2 } 



h- 



12 



where the scalar p = \p\, the spatial z = 1 ... 3, and the functions 

h(z ± )K (z ± ) - Ipjz^K^) hMKiM ^ 



are evaluated at z± = \p\jp\ + (p ± A 4 ) 2 - With these expressions it is explicitly verified 
that the normalization condition holds for any fi: 

1 = / (2^ ^ /i) = / (2^ K(~/^4(/i) + ^(-/i) ■ ^)] • (33) 
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